Mon. Not. R. Astron. Soc. 000, 000-000 (0000) 



Printed 1 February 2008 



(MN MfeK style file vl.4) 



Testing Tidal- Torque Theory: II. Alignment of Inertia and 
Shear and the Characteristics of Proto-haloes 



Cristiano Porciani 1 ' 2 , Avishai Dekel 1 and Yehuda Hoffman 

1 Racah Institute of Physics, The Hebrew University, Jerusalem 91904, Israel 

2 Institute of Astronomy, University of Cambridge, Madingley Road, Cambridge CB3-OHA, UK 



> 

in 



CN ! 1 February 2008 
O 
O 

CnI , ABSTRACT 

S3 , 

We investigate the cross-talk between the two key components of tidal-torque 
theory, the inertia (I) and shear (T) tensors, using a cosmological TV-body simulation 
C\) ' with thousands of well-resolved haloes. We find that the principal axes of / and T 

are strongly aligned, even though / characterizes the proto-halo locally while T is 
determined by the large-scale structure. Thus, the resultant galactic spin, which plays 
a key role in galaxy formation, is only a residual due to ~10 per cent deviations from 
perfect alignment of T and I, The T-I correlation induces a weak tendency for the 
proto-halo spin to be perpendicular to the major axes of T and /, but this correlation 
is erased by non-linear effects at late times, making the observed spins poor indicators 
ly-j . of the initial shear field. 

' However, the T-I correlation implies that the shear tensor can be used for identi- 

fying the positions and boundaries of proto-haloes in cosmological initial conditions — 
a missing piece in galaxy formation theory. The typical configuration is of a prolate 
proto-halo lying perpendicular to a large-scale high-density ridge, with the surround- 
ing voids inducing compression along the major and intermediate inertia axes of the 
proto-halo. This leads to a transient sub-halo filament along the large-scale ridge, 
whose sub-clumps then flow along the filament and merge into the final halo. 

The centres of proto-haloes tend to lie in ~ la over-density regions, but their 
association with linear density maxima smoothed on galactic scales is vague: only ~ 40 
per cent of the proto-haloes contain peaks within them. Several other characteristics 
distinguish proto-haloes from density peaks, e.g., they tend to compress along two 
principal axes while many peaks compress along three axes. 
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1 INTRODUCTION 

The angular-momentum properties of a galaxy system must 
have had a crucial role in determining its evolution and final 
type. A long history of research, involving Hoyle (1953) and 
Peebles (1969), led to a 'standard' theory for the origin of 
angular momentum in the cosmological framework of hierar- 
chical structure formation, the tidal-torque theory (hereafter 
TTT), as put together by Doroshkevich (1970) and White 
(1984). Special interest in the subject has been revived re- 
cently because of an 'angular-momentum crisis', arising from 
cosmological simulations of galaxy formation which seem to 
yield luminous galaxies that are significantly smaller, and 
of much less angular momentum than observed disc galax- 
ies (Navarro, Frenk & White 1995; Navarro & Steinmetz 
1997, 2000). Another current motivation comes from weak 



lensing studies, where the intrinsic distribution and align- 
ment of galaxy shapes, which may be derived from TTT, 
play an important role in interpreting the signal (Croft & 
Metzler 2000; Heavens, Refregier & Heymans 2000; Catelan, 
Kamionkowski & Blandford 2001; Catelan & Porciani 2001; 
Crittenden et al. 2001). 

These add a timely aspect to the motivation for re- 
visiting the classical problem of angular momentum, in an 
attempt to sharpen and deepen our understanding of its var- 
ious components and how they work together. In particular, 
this effort should start from the acquisition of angular mo- 
mentum by dark-matter haloes, despite the misleading ap- 
parent impression that this part of the theory is fairly well 
understood. The basic notion of TTT is that most of the an- 
gular momentum is being gained gradually by proto-haloes 
in the linear regime of density fluctuations growth, due to 
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tidal torques from neighboring fluctuations. This process is 
expected to continue as long as the proto-halo is expand- 
ing. Once it decouples from the expanding background and 
turns around to non-linear collapse and virialization, only 
little angular momentum is expected to be tidally exchanged 
between haloes. It is commonly assumed that the baryonic 
material, which in general follows the dark-matter distribu- 
tion inside each proto-halo, gains a similar specific angular 
momentum and carries it along when it contracts to form 
a luminous galaxy at the halo centre. This should allow us 
to predict galactic spins using the approximate but powerful 
analytic tools of quasi-linear theory of gravitational instabil- 
ity. Even if angular momentum is transferred from the gas to 
the dark matter during disc formation, the initial set-up of 
angular-momentum distribution in the halo is an important 
ingredient in the process. 

In a series of papers, we evaluate the performance of 
the TTT approximation, and trace the roles of its various 
ingredients and the cross-talk between them. We do it using 
a cosmological iV-body simulation with ~ 7300 well-resolved 
haloes. We find that some of the basic ingredients of stan- 
dard TTT involve certain unjustified assumptions and that 
our understanding of the theory is not full. These papers rep- 
resent attempts to clarify some of these controversial issues. 
In Paper I (Porciani, Dekel & Hoffman 2001), we evaluate 
how well does the approximation predict the final spin of a 
halo, given full knowledge of the corresponding initial proto- 
halo and the cosmological realization. In the present paper 
(Paper II), we attempt a deeper level of understanding of 
the origin of halo angular momentum, by investigating the 
relation between the different components of TTT. We find 
that this study connects to the fundamental open question 
of how to identify a proto-halo in a given realization of ini- 
tial conditions, and how proto-haloes evolve via collapse and 
mergers to the present virialized haloes. While the current 
discussion of proto-haloes and their evolution is limited to 
first results, a more detailed analysis will be reported in fu- 
ture papers. In Paper III (Porciani & Dekel in preparation), 
the standard scaling relation of TTT is revised, based on our 
finding that the shear tensor is weakly correlated with the 
density at the proto-halo centre of mass. This scaling rela- 
tion is used to predict the typical angular-momentum profile 
of haloes (Dekel et al. 2001; Bullock et al. 2001b), and to 
provide a simple way to incorporate angular momentum in 
semi-analytic models of galaxy formation (Mailer, Dekel & 
Somerville 2002). 

According to TTT, the angular momentum of a halo is 
due to the cross-talk between two key players: the inertia 
tensor /, describing the quadrupole structure of the proto- 
halo, and the shear tensor T, representing the external tidal 
field exerting the torque. There is a controversy in the lit- 
erature regarding the correlation between these two com- 
ponents, which makes a big difference in the outcome. In 
applications of TTT, it has commonly been assumed that 
I and T are uncorrelated (e.g. Hoffman 1986a,b; Heavens 
& Peacock 1988; Steinmetz & Bartelmann 1995; Catelan & 
Theuns 1996), based on the argument that the former is lo- 
cal and the latter must be dominated by external sources. 
Lee & Pen (2000, hereafter LP00), based on simulations 
with limited resolution, have raised some doubts concerning 
this assumption. Other applications have made the assump- 
tion that the present-day halo angular momentum correlates 



with the shear tensor of the initial conditions (LP00; Pen 
et al. 2000; Crittenden et al. 2001). This involves assum- 
ing both a strong correlation at the initial conditions and 
that this correlation survives the non-linear evolution at late 
times. In this paper, we investigate the validity of these as- 
sumptions, and find surprising results that shed new light 
on the basic understanding of TTT. 

One of our pleasant surprises is that a high degree of 
correlation between T and / leads to progress in a more 
general problem, that of identifying proto-haloes in the ini- 
tial conditions. Despite the impressive progress made in 
analysing the statistics of peaks in Gaussian random density 
fields (Bardeen et al. 1986; Hoffman 1986b; Bond & Myers 
1996), given a realization of such initial conditions we do 
not have a successful recipe for identifying the proto-halo 
centers and the Lagrangian region about them which will 
end up in the final virialized halo. We find that proto-haloes 
are only vaguely associated with density peaks. Straightfor- 
ward ideas for identifying the boundaries, such as involving 
iso-density or iso-potential contours (Bardeen et al. 1986; 
Hoffman 1986b, 1988b; Heavens & Peacock 1988; Catelan & 
Theuns 1996), do not provide a successful algorithm (van de 
Weygaert & Babul 1994; Bond & Myers 1996). These are key 
missing ingredients in a full theory for the origin of angu- 
lar momentum as well as for more general aspects of galaxy 
formation theory. We make here the first steps towards an 
algorithm that may provide the missing pieces. 

The outline of this paper is as follows: In §2| we summa- 
rize the basics of linear tidal-torque theory. In §3 we describe 
the simulation and halo finder, and the implementation of 
TTT to proto-haloes. In §[| we address the correlation be- 
tween the shear field and the inertia tensor of the proto-halo. 
In §JE]we show examples of the T — I correlation, and address 
the implications on the characteristics of proto-halo regions 
and how they evolve. In §^ we refer to other properties of 
proto-haloes. In §^ we address the correlation between the 
shear field and the spin direction. In §^| we discuss our results 
and conclude. 



2 TIDAL-TORQUE THEORY 

Here is a brief summary of the relevant basics of tidal-torque 
theory, which is described in some more detail in Paper I. 

The framework is the standard FRW cosmology in the 
matter-dominated era, with small density fluctuations that 
grow by gravitational instability. Given a proto-halo, a patch 
of matter occupying an Eulerian volume 7 that is destined 
to end up in a virialized halo, the goal is to compute its 
angular momentum about the centre of mass, to the low- 
est non- vanishing order in perturbation theory. The angular 
momentum at time t is 

L(t)= J p(T,t)[T(t)-T m (t)] X [v(i)-v cm (t)]d 3 r, (1) 

where r and v are the position and peculiar velocity, and 
the subscript cm denotes centre-of-mass quantities. Then, 
in comoving units x = r/a(t), 



L(t) = p(t)a 5 (t) / [l + {(x,t)][x(i)-x cm (f)]xx(()(fa!, 
J ~, 



(2) 
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where 5(x, t) is the density fluctuation field relative to the 
average density p(t), a(t) is the universal expansion factor, 
and the term proportional to v cm vanished. A dot denotes 
a derivative with respect to cosmic time t. 

The comoving Eulerian position of each fluid element is 
given by its initial, Lagrangian position q plus a displace- 
ment: x = q + S(q, t). When fluctuations are sufficiently 
small, or when the flow is properly smoothed, the mapping 
q — > x is reversible such that the flow is laminar. Then the 
Jacobian determinant J = ||9x/9q|| does not vanish, and 
the continuity equation implies f + <5[x(q, t)\ = J _1 (q, t). 
Substituting in Eq. (^|) one obtains 



L(t) = a 2 (t)p al 



[q-q + S(q,t)-S] X S(q,t)d 3 g, 



(3) 



where barred quantities are averages over q in V, the 
Lagrangian region corresponding to 7. The displacement 
S is now spelled out using the Zel'dovich approximation 
(Zel'dovich 1970), S(q,t) = -£>(t)V$(q), where $(q) = 
cf>(ci,t)/4nGp(t)a' 2 (t)D(t) (with G Newton's gravitational 
constant), and <j)(a,t) is the gravitational potential. Sub- 
stituting in Eq. (H) one obtains 



L(t) = -a\t)D(t) p al J (q - q) 



V$(q)d 3 g. (4) 



The explicit growth rate is L oc a 2 (t)D(t), which is oc t as 
long as the universe is Einstein-de Sitter, or close to flat 
(and matter dominated). 

Next, assume that the potential is varying smoothly 
within the volume F, such that it can be approximated br 
its second-order Taylor expansion about the centre of mass 



$(q') ~ $(0) + 



9$ 



1 d 2 $ 



2 dq' t dq' 



! f 

li 1j 



(5) 



where q' = q q. Substituting in Eq. (Q), one obtains the 
basic TTT expression for the ith Cartesian component: 



Li(t) = a (t)D(t) tijk T>ji lit ■ 



(6) 



where e^fc is the antisymmetric tensor, and the two key 
quantities are the deformation tensor at q' = 0, 



9 2 <E> 



dq'idq'j 



q' = 



and the inertia tensor of T, 



r - 3 
Hj = poa 



1 1 j3 / 



(7) 



(8) 



Note that only the traceless parts of the two tensors 
matter for the cross product in Eq. (^|). These are the ve- 
locity shear or tidal tensor, Tij = T>ij — (T>u/3)8ij, and the 
traceless quadrupolar inertia tensor, Iij — (Iu/3)Sij . Thus, to 
the first non-vanishing order, angular momentum is trans- 
ferred to the proto-halo by the gravitational coupling of the 
quadrupole moment of its mass distribution with the tidal 
field exerted by neighboring density fluctuations. The torque 
depends on the misalignment between the two. 



This is equivalent to assuming that the velocity field is well 



described by its linear Taylor expansion, i.e., u; 



3 TTT IN SIMULATIONS 

We summarize here the relevant issues concerning the sim- 
ulation, the halo finding, and the way we implement TTT. 
A more detailed description is provided in Paper I. 

The iV-body simulation was performed as part of the 
GIF project (e.g., Kauffmann et al. 1999) using the adap- 
tive P 3 M code developed by the Virgo consortium (Pearce et 
al. 1995; Pearce & Couchman 1997). As an example, we use 
a simulation of the rCDM scenario, in which the cosmology 
is flat, with density parameter f2 m — 1 and Hubble constant 
h — 0.5 (Ho — 100 h kms -1 Mpc -1 ). The power spectrum 
of initial density fluctuations is CDM with shape parame- 
ter r = 0.21, normalized to as = 0.51 today. The simulation 
was performed in a periodic cubic box of side 84.55 /i _1 Mpc, 
with 256 3 particles of 1.0 x lO lo /i _1 M each. Long-range 
gravitational forces were computed on a 512 3 mesh, while 
short-range interactions were calculated as in Efstathiou & 
Eastwood (1981). At late times (2 3), the gravitational 
potential asymptotically matches a Plummer law with soft- 
ening e = 36ft _1 kpc. The simulation started at z = 50 and 
ended at z = 0. The initial conditions were generated by dis- 
placing particles according to the Zel'dovich approximation 
from an initial stable 'glass' state (e.g. White 1996). More 
details are in Jenkins et al. (1998). 

Fig. [l] shows a typical halo and its proto-halo, to help 
illustrating the objects of our analysis. We identify dark- 
matter haloes at z = using a standard friends-of-friends 
algorithm with a linking length 0.2 in units of the aver- 
age inter-particle distance. This algorithm identifies regions 
bounded by a surface of approximately constant density, cor- 
responding to haloes with a mean density contrast ~ 180, in 
general agreement with spherical perturbations whose outer 
shells have collapsed recently. We then remove unbound par- 
ticles from each halo, and consider only haloes which contain 
more than 100 bound particles. Our results turn out to be 
insensitive to the removal of unbound particles, which is at 
a typical level of only a few per cent. 

The robustness of our conclusions regarding TTT with 
respect to the halo-finding algorithm (e.g., friends-of-friends 
versus fitting a spherical or ellipsoidal density profile, Bul- 
lock et al. 2001a) is investigated in another paper in prepa- 
ration. Note that all the haloes in our current sample are not 
subclumps of larger host haloes, for which linear theory is 
not expected to be valid. This excludes about 10 per cent of 
the haloes that are more massive than 10 12 h~ 1 M Q (Sigad 
et al. 2001). 

The proto-halo regions V are defined by simply tracing 
all the virialized halo particles, as identified today, into their 
Lagrangian positions. In most cases, most of the proto-halo 
mass is contained within a simply connected Lagrangian re- 
gion, but in some cases (~ 15 per cent) the proto-halo may 
be divided into two or three compact regions which are con- 
nected by thin filaments. Nearly 10 per cent of the proto- 
haloes are characterized by extended filaments departing 
from a compact core. A typical case is illustrated in Fig. |l|, 
where we show a halo at 2 = and its proto-halo at 2 = 50, 
embedded in the surrounding particles in a comoving slice 
that includes most of the halo particles. 

For each proto-halo, we compute the Lagrangian inertia 
tensor by direct summation over its N particles of mass m 
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Figure 1. An example of a halo at z = (center and right) and its proto-halo at z = 50 (left), embedded in the surrounding particles 
in a comoving slice about the halo centre of mass. The thickness of the slice is 5.4 and 0.9 1 Mpc (comoving) for the left, center and 
right panels respectively, such that all the halo particles are included. Particles belonging to the halo are marked by filled circles, while 
points denote background particles. Note that the halo is centred on one clump within a filament, while the neighbouring haloes are not 
marked. Note also that the proto-halo boundaries do not necesarily follow the orientation of the filament. Initial density contours in the 
region surrounding the proto-halo are shown in Figure 4. 



each: 

JV 

/^m^ftW^, (9) 

n=l 

with q the position of the n-th particle with respect to 
the halo centre of mass. 

The shear tensor at the proto-halo centre of mass is 
computed by first smoothing the potential used to generate 
the initial Zel'dovich displacements, and then differentiating 
it twice with respect to the spatial coordinates (Method 1). 
Smoothing is done using a top-hat window function, while 
derivatives are computed on a grid. The top-hat smoothing 
radius, in comoving units, is taken to be defined by the halo 
mass via (47r/3)po^ 3 = M. In Paper I we also tested two 
alternative methods for computing the shear tensor: Method 
2, in which the smoothing has been replaced by minimal 
variance fitting, and Method 3, also with minimal variance 
fitting, but in which we considered only the shear generated 
by the density perturbations lying outside the proto-halo 
volume. The TTT predictions for the halo spin based on the 
three methods were found to be of similar quality at z = 0, 
while Methods 2 and 3 are slightly more accurate at very 
high redshift. We adopt method 1 here because it does not 
depend explicitly on the detailed shape of the proto-haloes, 
which prevents spurious correlations between T and / due 
to the way T is computed. Also, this is the only method 
applicable in analytic and semi-analytic modelling. 

4 ALIGNMENT OF INERTIA AND SHEAR 

In Paper I, we evaluated the success of TTT in predict- 
ing the final halo angular momentum for a given proto-halo 



embedded in a given tidal environment at the initial con- 
ditions. For a deeper understanding of how TTT works, we 
now proceed to investigate the key players: the inertia and 
shear tensors and in particular the cross-talk between them. 
Note, in Eq. (^j), that angular momentum is generated only 
due to misalignments of the principal axes of these tensors. 
It has been assumed, in many occasions, that these two ten- 
sors are largely uncorrelated, the inertia tensor being a local 
property while the shear tensor is dominated by large-scale 
structure external to the proto-halo (e.g. Hoffman 1986a, 
1988a; Heavens & Peacock 1988; Steinmetz & Bartelmann 
1995; Catelan & Theuns 1996). Such a lack of correlation 
would have led to relatively large spins and would have pro- 
vided a specific statistical framework for TTT. On the other 
hand, a strong correlation between / and T would have led 
to relatively small spins, and would have invalidated some 
of the predictions based on the assumption of independence 
of I and T. Moreover, a correlation of this sort could pro- 
vide the missing clue for the special characteristics of the 
Lagrangian proto-halo regions F, those that will eventually 
evolve into virialized haloes. 

We address the correlation between the principal axes 
of / and T at the proto-halo centres of mass by first showing 
in Fig. [| the probability distributions of the cosines of the 
angles between them. The eigenvectors ij and t& are labeled 
in such a way that the corresponding eigenvalues are ranked 
ii > i 2 > is and ti < ti < £3, namely the major axes are 
denoted by 1 and the minor axes by 3 (and note, for ex- 
ample, that ti is the direction of maximum compression)]^] 
The highest spike (near a cosine of unity) indicates a very 

T The principal axes of Tij coincide with those of £>;j, because 
the difference of the two tensors is a scalar matrix. 
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Figure 2. Alignments of / and T. Distributions of the cosine of the angle between the directions of the principal axes of the inertia 
tensor (ij) and of the velocity shear tensor (tj.). The corresponding eigenvalues are ordered i\ > 12 > i'i and t\ < tt < tz, such that 1 
denotes the major axis and 3 denotes the minor axis. The first two moments of the distributions are given in Eq."Q and Eq. (jTlj). 



strong alignment between the minor axes of the two tensors, 
and the second-highest spike refers to a strong alignment 
between the major axes. The alignment between the inter- 
mediate axes is also apparent but somewhat weaker. The 
peaks near a cosine of zero indicate a significant tendency 
for orthogonality between the major axis of one tensor and 
the minor axis of the other. There is a somewhat weaker or- 
thogonality between the intermediate and minor axes, and 
the weakest orthogonality is between the major and inter- 
mediate axes. 

The first two moments of these distributions are 



IMj = • * j I > = 



and 



vu = <(ii - ij) > 



'0.838 0.363 0.106N 
0.370 0.793 0.196 
v 0.091 0.206 0.935/ 



'0.759 0.219 0.023^ 
0.225 0.692 0.083 
v 0.017 0.089 0.895/ 



(10) 



(11) 



Based on the higher-order moments and the number of 
haloes in the sample, the statistical uncertainty in these 
mean quantities is between 0.001 and 0.003, depending on 
the matrix element considered. In order to estimate the 
systematic errors due to the finite number of particles in 
each halo, we re-measured the moments considering only the 



haloes containing 1000 particles or more. The results agree 
with those dominated by smaller haloes of > 100 particles 
at the level of a few per cent. For example, we obtain 



/0.792 
Hij = 0.433 
1 0.117 



0.418 
0.733 
0.243 



O.U2\ 
0.226 
0.907 ; 



(12) 



with a typical statistical uncertainty of 0.010. 

The values of the matrix elements /xy and i/y can be 
used to quantify the degree of correlation between the shear 
and inertia tensors. Note that a perfect correlation corre- 
sponds to a diagonal, unit matrix both for fi and v, while 
perfect independence corresponds to all the elements being 
equal at 1/2 (for fj,) or 1/3 (for v). Note that, by definition, 
52 Vi i = ^Zi Ui i = 1- To characterize the above matrices 
with single numbers, we define a correlation parameter for 
fiij by 



1 



(Mi ~ I/ 2 
-1/2 



- 1 2 

- 3 + g ijMi 



(13) 



it obtains the values zero and unity for minimum and max- 
imum correlation respectively. Similarly, we define a corre- 
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Figure 3. Examples of the correlation between I and T. The top panels show maps of the density (left) and velocity (right) fields, at 
z = 50, in a section of an X — Y plane from the simulation box. The fields are smoothed with a top-hat window of radius 0.95ft — 1 Mpc 
corresponding to a 100-particlc halo. The contours refer to the density contrast linearly extrapolated to z = 0. The bottom panels show 
all the proto-haloes whose centres of mass lie within one smoothing length of the plane (A cluster-size halo associated with the high 
density peaks at the top is not shown because its centre lies further away from the plane). The left panel shows the projection of the 
proto-halo particle positions. The haloes contain, from left to right, 179, 257, 143, 168, 100, 183, and 300 particles. The right panel 
shows the projections of the major axes of / (dark lines) and T (light lines) about the centres of mass (filled circles). The line length is 
proportional to the projection of a unit vector along the major axis on to the X — Y plane. To set the scale, we note, for example, that 
the principal axis of the inertia tensor of the proto-halo at the centre of the panel lies almost exactly in the plane shown. 



lation parameter using Vij, 

ly gj ~ 1/3 = 1 ( _ 1} (14) 
9 ^ Sij - 1/3 2 v ; y ' 

We find in the simulation c M = 0.61 and c„ = 0.67, indicating 
a strong correlation between I and T. 

We have repeated the above correlation analysis for the 
two alternative methods of computing the shear tensor de- 
scribed in §4 of Paper I. When we minimize smoothing and 
include only the part of the shear tensor that is generated 
by fluctuations external to the proto-halo (Method 3), the 
correlations are very similar to those obtained with top-hat 



smoothing, indicating that these are indeed the correlations 
of physical significance that we are after. When we minimize 
smoothing and include the contribution from internal fluc- 
tuations (Method 2), the correlations still show but they get 
significantly weaker because of the additional noise in the 
computation of T. Thus, our standard way of computing T 
(Method 1) seems to pick up properly the correlation with 
I, and we proceed with it. 

Note that T-I correlations are unavoidable in the frame- 
work of the Zel'dovich approximation, where a halo consists 
by the matter that has shell-crossed. However, the strength 
of the expected correlations has still to be worked out. 
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The implication of the T-I alignment is that proto-halo 
spins are due to small residuals from this correlation. This 
means that the haloes acquire significantly less angular mo- 
mentum than one would have expected based on a simple 
dimensional analysis that ignores the correlation. We quan- 
tify the effect by a comparison to an artificial case where 
the alignments are erased. For each halo in our sample, we 
randomize the relative directions of the eigenvectors of the 
inertia and tidal tensors (by a three-dimensional rotation 
with random Euler angles), and then re-compute the angu- 
lar momentum using Eq. (p|). A vast majority of the haloes, 
about 88 per cent, are found to be associated with a larger 
spin after the randomization procedure. The T-I correla- 
tion is found to reduce the halo spin amplitude with respect 
to the randomized case by an average factor of 3.f, with 
a comparable scatter of 3.2. This heads towards explaining 
why haloes have such low values of the dimensionless spin 
parameter, A ~ 0.035 on average (Barnes & Efstathiou f987; 



Bullock et al. 2001b), namely a very small rotational energy 
compared to gravitational or kinetic energy. 



5 PROTO-HALO REGIONS 

The strong correlation between the principal axes of the in- 
ertia and tidal tensors promises very interesting implications 
on galaxy-formation theory. It indicates that the tidal field 
plays a key role in determining the locations and shapes of 
proto-haloes, and therefore may provide a useful tool for 
identifying proto-haloes in cosmological initial conditions. 

As a first step, we try to gain intuitive understanding of 
the nature of the correlation between T and I by inspecting 
a few proto-haloes and their cross-talk with their cosmolog- 
ical environment. Fig. ^ and Fig. |i| show a few examples, 
in regions of the simulation box containing several proto- 
haloes. In the top panels we show maps of the density and 
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velocity fields at the initial conditions, to indicate some of 
the qualitative properties of the shear field. We focus on a 
section of one plane at at time, and project the velocity field 
onto it. The fields are smoothed with a top-hat window cor- 
responding to the typical halo masses shown in that region: 
in Fig. ^, the smoothing radius is 0.95 /i _1 Mpc, to match the 
haloes that contain 100 to 300 particles each, while in Fig. ^ 
it is 2.58 /i _1 Mpc, to match the haloes of 1000 to 2000 parti- 
cles. The smoothed density contrast is linearly extrapolated 
to z = 0. 

We then introduce, in the bottom panels, the proto- 
haloes and their inertia tensors. We present all the proto- 
haloes whose centres of mass lie within one smoothing length 
of the plane by showing, in the bottom-left panels, the pro- 
jection of the proto-halo particle positions onto the plane. 
We see a tendency for an association between the proto-halo 
centres of mass and density peaks smoothed on a halo scale, 
but this association seems to be far from perfect (e.g. near 
the centre of the frame in Fig. |^) . Note that the single plane 
shown does not allow an accurate identification of the peak 
location, and therefore a definite evaluation of the associa- 
tion of proto-haloes and peaks cannot be properly addressed 
by observing these maps. In order to quantify the correlation 
between proto-haloes and peaks, we identified peaks in the 
linear density field, smoothed with a top-hat window that 
contains n (with n=100, 1000 and 10000) particles, and de- 
termined the nearest peak to the centre of mass of each of 
the proto-haloes that contain n (±10 per cent) particles. In- 
dependently of the halo mass, we find that only ~ 35 — 45 
per cent of the proto-halo centers lie within one smooth- 
ing radius from the nearest peak, and ~ 60 — 65 per cent 
lie within two smoothing radii. We conclude that while the 
proto-haloes tend to lie in high density regions, their centres 
do not coincide with the local density maxima, and in fact 
the spatial correlation between the two is quite weak. 

In the bottom-right panels, we stress the projections 
of the major axes of the inertia and shear tensors about the 
centres of mass, and indicate by the length of the line the co- 
sine of the angle between the axis and the plane shown. Note, 
for example, that the major axes of those haloes in Fig. ^ 
that reside near the centre of the frame and towards the bot- 
tom right happen to lie almost perfectly in the plane shown. 
All the proto-haloes that show a significant deviation from 
spherical symmetry have their major inertia axis strongly 
aligned with the first principal shear axis at their centre of 
mass. This alignment reflects the fact that the largest com- 
pression flow towards the centre of the proto-halo is along 
the major inertia axis of the proto-halo. Indeed, this is ex- 
actly what is required in order to compress the elongated 
proto-halo into the more centrally concentrated and quite 
spherical configuration identified by the halo finder after col- 
lapse and virialization. 

The detected alignment is a clear manifestation of the 
crucial role played by the external shear in determining the 
shape of the Lagrangian volume of the proto-halo. One way 
to interpret this alignment is as follows. The compression 
along the major axis of T is associated with a flow of matter 
from the vicinity into the proto-halo. This matter comes 
from relatively large distances and it is therefore responsible 
for a large inertia moment along this axis in the proto-halo 
configuration. In contrast, the dilation along the minor axis 
of T causes matter to be tidally stripped, thus leading to a 



relatively small inertia moment of the proto-halo along this 
axis. This picture, in which the boundaries of the proto- 
halo are fixed by the push and pull of the external mass 
distribution, is in some sense the opposite of the common 
wisdom, where the crucial factor is assumed to be the self- 
gravity attraction. This other approach, which predicts that 
collapse first occurs along the minor inertia axis, is based on 
the simplified model of the collapse of an isolated ellipsoidal 
perturbation that starts at rest or comoving with the Hubble 
flow (Lin, Mestel & Shu 1965; Zel'dovich 1965). 

Snapshots of the evolution of a typical halo (the same 
massive halo shown in Fig. |l| and Fig. are shown in Fig. ^[ 
Plotted are the positions at different epochs of the particles 
that form the halo identified at z = 0. The proto-halo first 
collapses along its major and intermediate inertia axes, giv- 
ing rise to an elongated structure made of sub-clumps. The 
final halo is then assembled by merging and accretion along 
the axis of the elongated filament. This transient sub-halo 
filament lies along the large-scale filament in which the fi- 
nal halo is embedded (see Fig. [I] and Fig. ^J). It indicates 
that the late stages of halo formation are associated with 
flows along preferential directions aligned with the cosmic 
web (see also Colberg et al. 1999). 

It has often been assumed that the boundaries of proto- 
haloes could be identified with some threshold iso-density 
contours about density maxima. (Bardeen et al. 1986; Hoff- 
man 1986b, 1988b; Heavens & Peacock 1988; Catelan & 
Theuns 1996). In this case, the density distribution in a 
proto-halo could be approximated by a second-order Taylor 
expansion of the density profile about the peak. Figures ^ 
and ^ demonstrate that this assumption is going to fail in 
most cases shown (except, perhaps, the halo at the bottom- 
right of Fig. ^). In a typical proto-halo, the boundaries are 
not determined by the self gravity due to the local mass dis- 
tribution within the proto-halo but rather by the tidal field 
due to the external mass distribution. Such a proto-halo is 
typically embedded within a large-scale elongated density 
ridge, but is not necessarily centred on a local density peak. 
The compressions exerted by the void regions surrounding 
the ridge define the compression axes of the shear tensor in 
directions perpendicular to the ridge, while the tides from 
the other parts of the ridge cause large-scale dilation along 
the ridge. The associated push and pull of mass along these 
directions make the large inertia axes of the proto-halo lie 
perpendicular to the ridge, and the minor inertia axis lie 
parallel to the ridge. This understanding should be trans- 
lated to a practical recipe for identifying the boundaries of 
proto-haloes. At later stages we also see an internal com- 
pression and merging along the filament, namely the minor 
inertia axis of the proto-halo. 



6 OTHER PROPERTIES OF PROTO-HALOES 

We continue with a complementary investigation of the 
properties of proto-haloes, using several different statistics 
that may, in particular, distinguish them from random Gaus- 
sian peaks. 

First of all, we study the shapes of the proto-haloes 
in Lagrangian space using their inertia tensors. A complete 
description of the statistical properties of the eigenvalues 
of the inertia tensor (a symmetric, positive definite matrix) 
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Figure 5. Time evolution of a protohalo. Shown are the projected positions of the halo particles at different redshifts in a fixed comoving 
box about the centre of mass. The first collapse along the major axis of inertia (and the major axis of the tidal field) leads to the formation 
of a transient filament, which breaks into sub-clumps. These sub-clumps then flow along the filament and merge into the final halo 



can be obtained in terms of the three following parameters: 
the trace r = ii + «2 + is, the ellipticity e = (ii — is)/2r, 
and the prolateness p = (i% — 2i2 + iz)/2r. A perfect sphere 
has e = p = 0. A flat circular disc (ultimate oblateness) 
has e = 1/4 and p = —1/4. A thin filament (ultimate pro- 
lateness) has e — 1/2 and p = 1/2. Thus, e measures the 
deviation from sphericity and p measures the prolateness 
versus oblateness. In Fig. |we show the joint distribution of 
ellipticity and prolateness, and the probability distribution 
marginalized along each axis, for our proto-haloes at z — 50. 
The boundaries e > —p, e > p, and p > 3e — 1 arise from 
the conditions ii > 12, «2 > 13, and > 0, respectively. As a 
consequence, the data-points populate only a triangle in the 
e-p plane, with vertices at (0, 0), (1/2, 1/2) and (1/4, -1/4). 
This introduces a correlation between e and p at high ellip- 
ticities (e > 1/4). We find that only 3 per cent of the proto- 
haloes are nearly spherical, with e < 0.1 (or \p\ < 0.1). Most 
of the proto-haloes, 74 per cent, have moderate but signif- 
icant ellipticities in the range 0.1 < e < 0.25. A significant 
fraction, 23 per cent, have extreme ellipticities of e > 0.25 
(and almost all are prolate configurations). About 68 per 
cent of the proto-haloes are prolate, p > 0, and about 2.3 per 
cent are extremely prolate, p > 0.25. The average and stan- 
dard deviation for e and p are 0.206±0.062 and 0.048±0.092, 
respectively. It is pretty surprising to note that, despite the 



pronounced difference of the proto-halo boundaries from iso- 
density contours about density peaks, the shape distribution 
of the two are not very different. Bardeen et al. (1986, sec- 
tion VII) quote for overdensity patches about 2.7a peaks: 
0.17 ± 0.07 and 0.005 ± 0.098 for e and p respectively. Using 
their equations 7.7 and 6.17 for la peaks, which we found 
to be more appropriate for proto-haloes, we get 0.19 ± 0.08 
and 0.04 ± 0.11 for e and p respectively. 

In §^ we found that the principal axes of the inertia and 
tidal (and deformation) tensors tend to be aligned, which 
means that proto-haloes will have their dominant collapse 
along their major axes of inertia. We now test whether 
the shape of the proto-halo is correlated with the relative 
strength of the eigenvalues of the deformation tensor, i.e., 
the compression factors along the principal axes of T>^ In 
Fig. 0, we show the joint distributions of the ratios ij/ik 
and dj/dk between the eigenvalues of the inertia and defor- 
mation tensors (dj = ti + X>ii/3)|j] Certain correlations are 

• We deal here with T> rather than T because we expect the self 
gravity of the proto-halo (described by the diagonal terms of T>) 
to also play a role in shaping it, as in the basic formulation of the 
Zel'dovich approximation. 

§ We prefer here these ratios over the e and p parameters used to 
characterize / because the latter are ill-defined for non-positive- 
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Figure 6. Prolateness and ellipticity for the inertia tensor of 
proto-haloes. The average and standard deviation values are (e) = 
0.206 and <r e = 0.062 and (p) = 0.048 and a v = 0.092. 
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detected for some of these ratios, meaning that the degree 
of asymmetry between the principal axes of the tidal field 
plays a role in determining the Lagrangian shape of a proto- 
halo. For instance, it appears that the ratio Gfo/di is anti- 
correlated with v$jii. This may be interpreted as another 
face of the T-I correlation (§^), where the largest (smallest) 
compression flow in a proto-halo is associated with its major 
(minor) axis of inertia. For example, when the two largest 
compression flows are of comparable strength, di ~ da > d,3, 
the proto-halo configuration is oblate, 13 < 12 ~ ii, namely 
a large c^/di and a small izjii- On the other hand, when 
the compression is mainly along one axis and c?2 — (i.e. 
small dz/di), the configuration is prolate, 13 ~ ii < ii (i.e. 
large 13/12). Note that these correlations, much like the T-I 
correlation, derive from the constraint that the Lagrangian 
region is destined to end up in a virialized halo. When pan- 
cakes or filaments form instead, the expected correlations 
may be different, with the final minor axis of inertia corre- 
lated with the direction of the largest flow of compression, 
etc. As seen in Fig. ^, the typical evolution from an initial 
proto-halo to a final halo indeed passes through an interme- 
diate phase of a pancake (or a filament) lying perpendicular 
to the direction(s) of first collapse. 

Back to proto-haloes, it is clear that binary correlations 
cannot tell the full story; for a complete study one should 
consider the joint distribution of the three eigenvalues. Other 
quantities, such as bulk flows, may also influence the shapes 
of proto-haloes. Therefore, this issue deserves a more de- 
tailed study beyond the scope of the current paper, which 
only provides first clues. 

From the dynamical point of view, proto-haloes can also 
be classified by the signs of the eigenvalues of the deforma- 
tion tensor at the centre of mass, indicating whether the ini- 
tial flows along the principal directions are of compression of 
dilation. We find that a small minority of the proto-haloes, 
about 11 per cent, are contracting along all three principal 
axes. The vast majority, about 86 per cent, are initially col- 
lapsing along two directions and expanding along the third. 
Only 2.7 per cent are collapsing along one direction and ex- 
panding along two.p| This clearly distinguishes the centres 
of mass of proto-haloes from random points in the Gaussian 
field, where the probability of the corresponding dynamical 
configurations based on the deformation tensor would have 
been 8, 42, and 42 per cent respectively, with the rest 8 per 
cent expanding along three directions (Doroshkevich 1970). 
It also distinguishes the proto-haloes from peaks of the lin- 



Figure 7. Ratios of the eigenvalues of the inertia tensor versus 
ratios of the eigenvalues of the deformation tensor (evaluated at 
the halo centre of mass) for the entire halo population in the 
simulation. Also shown is the corresponding linear correlation co- 
efficient for the joint distribution in each panel. 



definite matrices such as T>, for which the trace and the single 
eigenvalues may vanish or be negative. 

^ These fractions are obtained with standard top-hat smoothing 
of the deformation tensor. Similar numbers (15, 83, 1.5 per cent) 
are found when only the external velocity field is taken into ac- 
count and minimal smoothing is applied (Method 3 of Paper I). 
However, when the deformation tensor is computed with minimal 
smoothing and includes the contribution of fluctuations inside 
the proto-halo (Method 2), the frequencies become 98.3, 1.7, 0.0 
per cent, namely almost all the proto-haloes are collapsing along 
three spatial dimensions. This is because the local gravitational 
attraction towards the centre is dominant over the external shear 
such that it turns the large-scale expansion associated with d% 
(and 13) into a local contraction. This could be a feature that 
distinguishes proto-haloes from random patches. 
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Figure 8. Proto-haloes versus peaks: Distribution of density 
height v = <5/<t at proto-halo centres (solid histogram) versus 
random points (dotted line) and density maxima (dashed line). 
The density is smoothed top-hat over the proto-halo scale. The 
predicted probability distributions for density maxima and for 
random points are for a corresponding linear Gaussian overdensity 
field, smoothed on the scale of the smallest halo in the sample. 
The main panel refers to all the haloes containing more than 
100 particles, while the inset refers to haloes of more than 1000 
particles. The analytic estimates for peaks practically coincide 
with the actual density distribution of peaks in our simulation. 
The averages values, of 1.05 and 1.35, are marked by the dot- 
dashed lines. The corresponding standard deviations are 0.37 and 
0.33, respectively. 



ear density field. When smoothed with a top-hat window 
corresponding to 100 particles, we find in our simulation 
that the probabilities of the three dynamical configurations 
are 45, 46 and 8.6 per cent, respectively. These fractions 
somewhat depend on the smoothing length. Using a top-hat 
window containing 1000 (10000) particles, the corresponding 
probabilities become 50 (67), 44 (30) and 5.4 (3.4) per cent. 
The peaks are characterized by compression along three or 
two directions, while most of the proto-haloes are compress- 
ing along two directions. This finding is consistent with the 
picture arising from §|H] where proto-haloes are typically em- 
bedded in elongated, filament- like large-scale structures, sur- 
rounded by voids that induce compression flows along the 
two principal directions orthogonal to the filament. 

A certain fraction of our proto-haloes are associated, 
to some degree, with density peaks in the initial conditions. 
In Fig. ^ we show the distribution of v = 5 /a over all the 
haloes containing more than 100 (main panel) or 1000 (in- 
set) particles, where S is the density contrast at the centre 
of mass of the proto-halo and a is the rms density contrast 
over all space. This is compared with the expected distri- 
bution of v, for random points and for density peaks, in 
a Gaussian random field with the rCDM power spectrum, 
smoothed on the scale corresponding to 100 (main panel) 
and 1000 (inset) particles of our simulation. Note that our 
proto-haloes typically correspond to one-<r fluctuations with 



Figure 9. Comparison to the spherical collapse model: Distri- 
bution of overdensity for proto-haloes, evaluated at the centre of 
mass, smoothed top-hat on the halo scale, and linearly extrapo- 
lated from the initial conditions to z = 0. The average value of 
the distribution, (S) = 1.55, is marked by the dashed line. The 
standard deviation of the distribution is 0.54. The value predicted 
by the spherical collapse model, 8 ~ 1.686, is marked by the dot- 
ted line. For haloes with more than 1000 particles, the average is 
1.47 and the standard deviation is 0.36. 

a relatively small dispersion. This distribution can be clearly 
distinguished from the density distribution in randomly se- 
lected density maxima, which average at about 1.4(7, and 
have extended tails - especially towards high values. Some 
of the high peaks are embedded in more extended pertur- 
bations of v ~ 1. They give rise to clumps that eventually 
merge into the larger haloes that we identify today, and they 
are therefore not included in our sample of proto-haloes. 
Also, recall that the impression from Fig. ^ and Fig. [| was 
that the correlation between the halo centres of mass and 
the positions of density peaks is quite limited, and therefore 
the central proto-halo densities tend to be lower than the 
heights of the associated peaks. It is also possible that some 
very high and very low peaks may have deformation config- 
urations that would end up as very aspherical or fragmented 
structures rather than coherent virialized haloes (e.g. Katz, 
Quinn & Gelb 1993; van de Weygaert & Babul 1994). 

One often uses the terminology of the spherical col- 
lapse model to characterize the evolution of proto-haloes. 
For instance, one refers to 'turnaround' and 'collapse' time, 
even though, strictly speaking, these quantities are not well- 
defined for non-spherical objects (see, however, the discus- 
sion in Sugerman et al. 2000) . One way to test our definition 
of just-collapsed haloes at z = 0, and to address the spheric- 
ity of proto-haloes, is by evaluating the accuracy with which 
the spherical collapse model describes their evolution. Fig. ^ 
shows the distribution of the linearly extrapolated density 
contrast at the proto-halo centre of mass, smoothed on the 
halo size. The average of the distribution, S — 1.55, is within 
~10 per cent of the standard prediction of the spherical col- 
lapse model at collapse time, S ~ 1.686. This indicates that 
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Figure 10. Correlation between halo spin and the linear shear 
tensor. Shown is the distribution of the cosine of the angle be- 
tween the halo spin (at z = 50 or at z = 0) and each of the three 
eigenvectors of the linear shear tensor. The averages and standard 
deviations at z = 50, for axes 1,2,3 respectively, are: 0.35 ± 0.27, 
0.56 ± 0.30, 0.55 ± 0.31. The corresponding values at z = are: 
0.52 ± 0.30, 0.50 ± 0.30, 0.47 ± 0.29. 

our halo finding method and the spherical model are con- 
sistent on average. However, the scatter is large, indicating 
strong deviations from sphericity for many individual proto- 
haloes, which implies a big uncertainty in the turnaround or 
collapse times. 

What are the proto-haloes that populate the tails of the 
linear overdensity distribution? In particular, how do proto- 
haloes of 8 <C 1 manage to make haloes by today? It is ex- 
pected that the presence of shear may speed-up the collapse 
of a density perturbation with respect to the spherical case 
(Hoffman 1986b; Zaroubi & Hoffman 1993; Bertschinger & 
Jain 1994). In fact, our haloes with 8 <C 1 are indeed charac- 
terized by a strong shear. They typically have di ~ —d% and 

<C |d3 1, i.e. they are characterized by compression and 
dilation factors of similar amplitudes on orthogonal axes. 
In this case, shear terms in the Raychaudhuri equation (e.g. 
Bertschinger & Jain 1994) account for significant corrections 
to the spherical terms, already at z = 50. On the other hand, 
there are high density proto-haloes, S ~ 3, which collapse 
only today, which is late compared to the predictions of the 
spherical model. These are generally characterized by strong 
compression flows along 2 directions and by mild expansion 
along the third axis. The reason for their late collapse is not 
obvious and deserves further investigation which is beyond 
our scope here. 



7 ALIGNMENT OF SPIN AND SHEAR 

Had the inertia and shear tensors been uncorrelated, one 
could have argued that the direction of the spin in the lin- 
ear regime should tend to be aligned with the direction of 



the middle eigenvector of T (and of 7). The Cartesian com- 
ponents of L in the frame where T is diagonal are 

L l tx (tj - t k )Ijk , (15) 

where i,j,k are cyclic permutations of 1,2,3. We average 
over all the possible orthogonal matrices R that may relate 
the uncorrelated principal frames of T and I, while keeping 
the eigenvalues of the two tensors fixed, and obtain 

(\Li\) <x \tj - t k \ { \Ij k \ \i\,ii,i3) = \tj - tk\ f{ii,i2,h) ■ 

(16) 

The assumed independence of T and I would have guaran- 
teed that in the principal frame of T the conditional aver- 
age of Ijk (given the eigenvalues of /) would have been the 
same for all j 7^ k, namely a certain function f(ii,i2,is). 
On average, the largest component \Li\ would have been 
the one for which \tj — tk\ is the largest, which is necessarily 
L2 oc I £3 — t\\, because by definition ti < ti < £3. Thus, 
the angular momentum would have tended to be aligned 

with t2.P] To demonstrate this argument and evaluate the 
strength of the effect, we use the randomized T-I pairs ob- 
tained as described at the end of In this case, we find 
(L ■ ti) = 0.40,0.66,0.38 (for i = 1,2,3), with a scatter of 
~ 0.3 about these mean values. We see that the tendency for 
alignment with the middle eigenvector is quite significant. 

Given that the inertia and shear tensors are in fact 
strongly correlated (§fl), we now investigate, using TTT, 
what kind of alignment, if any, one should expect between 
L and T at the initial conditions. Due to the T — I corre- 
lation, the distribution of the orthogonal linear transforma- 
tions R should not be uniform but rather favour rotations 
with small Euler angles. Moreover, since we saw in that 
different pairs of eigenvectors have different degrees of align- 
ment, the average off-diagonal terms of I in the principal 
frame of T would not be all equal anymore. This could be 
sufficient for drastically changing the preferential alignment 
of the angular momentum with t2 seen in the uncorrelated 
case. Using TTT based on the inertia and tidal tensors for 
proto-haloes as extracted from the simulation, we actually 
obtain (L • tj) = 0.34,0.59,0.52, with a scatter of ~ 0.3 in 
every case. This shows that the T-I correlation indeed mod- 
ifies the L-T alignment predicted in the uncorrelated case. 
The spin tendency for alignment with t2 is still the strongest, 
but it is weaker than before. The alignment with t3 becomes 
stronger and almost comparable to the alignment with ta, 
while the alignment with ti becomes even weaker than it 
was, namely the spin prefers to be perpendicular to ti. We 
see that, generally speaking, the correlation between I and 
T tends to weaken the correlation between L and T. 

Next, we measure the actual spins of the proto-haloes in 
the initial conditions of the simulation, and correlate it with 
the directions of the eigenvectors of T. The distributions 
of angles between these vectors are shown in Fig. [l^ (top 
panels). We see that the spin has a weak but significant 
tendency to be perpendicular to the major axis of T, and an 
even weaker tendency to be parallel either to the middle or 

II Note that a similar argument can be used to show that L would 
have also been preferentially aligned with the middle eigenvector 
of I. In this case, however, since all the eigenvalues ij are positive, 
one expects a weaker correlation. 
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to the minor axis of T. In this case, (L-ti) = 0.35, 0.56, 0.55, 
in good agreement with the TTT approximation. 

Thus, contrary to LP00, we find no evidence for a strong 
correlation with the direction of the middle axis of T. LP00 
have proposed to parametrize L-T correlations through the 
expression 

(L i Li\f ij ) = ^5 ij -af ik f kj , (17) 

where Tij = Tij /[TikTki] 1 ^ 2 is the trace-free unit shear ten- 
sor, and a is the correlation parameter. This is particu- 
larly convenient for analytical studies of spin statistics (e.g. 
Crittenden et al. 2001), since it bypasses the ambiguous is- 
sue of determining the inertia tensors of proto-haloes and 
their statistics. The correlation parameter a vanishes if L is 
randomly distributed with respect to T, and according to 
LP00 it should take the value ~ 3/5 if TTT holds and T is 
independent of I (empirically, we actually get a — 0.68±0.01 
from the randomized sample in our simulation). It is espe- 
cially interesting to provide an independent measure of a 
from our simulation because Eq. ( p"7j ) has been used in mod- 
eling spin-spin correlations (Crittenden et al. 2001) and only 
one determination of a is available (LP00). 

The correlation parameter a can be evaluated in the 
principal frame of the tidal tensor. Denoting the eigenvalues 
and eigenvectors of T by t, and t<, one obtains a = 2 — 

6(Ei(L ■ t0 2 ^>. wnere E*i = V 2 has be en used (LP00). 
Using this expression for a, we obtain by averaging over all 
the proto-haloes in our initial conditions a — 0.28±0.01, and 
when we restrict ourselves to proto-haloes containing 1000 
particles or more we get a — 0.33±0.04. It is worth stressing, 
however, that the dispersion among the proto-haloes of the 
quantity that is averaged in the determination of a is very 
large, with an rms value of about 0.9, which is not a good 
property of this statistic. 

Next, we wish to find out how much of the (weak) 
L-T correlation detected at the initial conditions actually 
survives the non-linear evolution to the present epoch. For 
this we plot in the bottom panels of Fig. |l(j the distribu- 
tions of angles between the final spin and the initial shear 
tensor. We find that the correlation is strongly weakened. 
This is hardly surprising based on the significant non-linear 
evolution in spin direction found in Paper I. We obtain 
a — 0.07 ± 0.01, at variance with the result obtained by 
LP00, a = 0.24 ± 0.02 at z = 0, using a low-resolution 
particle-mesh simulation of an fl — 1 CDM model. Our re- 
sult persists when considering only the haloes containing 
1000 particles or more, a = 0.07 ± 0.04. The discrepancy 
with LP00 is probably an effect of the higher accuracy with 
which our high-resolution, adaptive P 3 M code accounts for 
non-linear effects. The smoothing length used to define Tij 
in the simulation could also affect the results (for instance, 
smoothing on too small a scale might give unphysically small 
L — T correlations). In order to test the dependence of our 

** We notice that the claim of LP00, that Eq. ( |l7| ) represents the 
most general quadratic relation between a unit vector and a unit 
tensor, is not formally true, because, for example, Eq. ( |l7[ ) cannot 
describe the case where L systematically lies along one of the 
principal axes of T. However, since this configuration cannot be 
realized in TTT, we adopt Eq. ([n]) as a plausible parametrization 
of the correlation. 



findings on the adopted filtering radius, we re-compute a 
using a top-hat window function which contains 8 times the 
halo mass. In this case, for haloes with 100 particles or more, 
we find a = 0.24±0.01 at z=50, and a = 0.08±0.01 at z = 0. 
This shows that smoothing is not a major issue here. 

Our result above seems to suggest that observed spin di- 
rections, as deduced from disc orientations, cannot be used 
to deduce useful information about the initial shear tensor. 
How relevant are our results to the orientation of disc galax- 
ies? On one hand, the haloes in our current analysis are all 
characterized by a present-day overdensity of ~ 180, while 
the galactic discs are associated with much higher density 
contrasts, so a certain unknown extrapolation is required. 
Moreover, it is plausible that the plane of today's disc galaxy 
is determined by the halo spin at an early time, before the 
gas fell into the disc. As expected, we see in the simulation 
that the non-linear effects wipe out the spin-shear correla- 
tion in a gradual fashion, e.g., a — 0.26,0.20,0.11 (±0.01) 
at z — 1.2,0.6,0.2 respectively. Therefore, the disc orienta- 
tion may after all preserve some of the weak correlation with 
the initial shear field, and this should be true especially for 
galaxies that formed at high redshift. On the other hand, 
we know that the spin and shear are correlated better at 
high redshift only for the haloes that we selected at z — 0, 
and we do not know whether this is true for the subhaloes 
which host disc formation at higher redshift. It is also pos- 
sible that the statistical orthogonality between the spin and 
the first principal shear axis is preserved better in specific re- 
gions were cosmic flows are particularly cold. The cold-flow 
neighborhood of the Local Group might provide a suitable 
field for studying such a correlation. Still, we expect only a 
weak correlation. 



8 CONCLUSION 

In order to deepen our understanding of how tidal torques 
actually work, we investigated the cross-talk between the 
main components of the theory using ~ 7300 well-resolved 
haloes extracted from a cosmological JV-body simulation. 
We studied the correlation between the proto-halo inertia 
tensor and the external shear tensor, and between them and 
the spin direction, and tried to characterize the proto-halo 
regions in several different ways. We defined haloes in to- 
day's density field using the standard friends-of-friends al- 
gorithm, but it would be useful for future work to check 
robustness to alternative halo finders as well as to different 
cosmological scenarios. 

We found to our surprise that the T and I tensors are 
strongly correlated, in the sense that their minor, major and 
middle principal axes tend to be aligned, in this order. This 
means that the angular momentum, which plays such a cru- 
cial role in the formation of disc galaxies, is only a residual 
which arises from the little, ~10 per cent deviations from 
perfect alignment of T and /. 

The T — I correlation induces a weak tendency of the 
proto-halo spin to be perpendicular to the major axis of 
T (and /). It also slightly weakens the spin tendency to 
be aligned with the middle axis (a tendency that would 
have dominated had T and I been uncorrelated) and slightly 
strengthen its alignment with the minor axis. However, these 
correlations, even at the initial conditions, are weak. Fur- 
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thermore, non-linear changes in spin direction at late times 
practically erase the memory of the initial shear tensor, and 
therefore observed spin directions cannot serve as very use- 
ful indicators for the initial shear tensor (again, in variance 
with LPOO). The only partial caveat is that today's discs 
may reflect the halo spin directions at some high redshift, 
which may still preserve some weak correlation with the lin- 
ear shear field. A study of this effect would require a high- 
resolution simulation in which a detailed galaxy formation 
scheme is incorporated. 

On the other hand, the strong T — I correlation provides 
a promising hint for how to make progress in a long-standing 
open question in galaxy formation theory. That is, how to 
characterize proto-haloes and their boundaries. Speculations 
based on iso-density or iso-potential contours about high- 
density peaks have failed the tests of simulations. We first 
realize that the centres of proto-haloes tend to lie in ~ la 
over-density regions, but their association with the linear 
density maxima smoothed on galactic scales is quite limited. 
We find that the proto-haloes tend to be elongated along the 
direction of maximum compression of the smoothed velocity 
field. A typical configuration is of an elongated proto-halo ly- 
ing perpendicular to an elongated background density ridge, 
with neighboring voids inducing the compression along the 
major and intermediate inertia axes. The T — I correlation 
detected here should thus enable the construction of a de- 
tailed algorithm to identify proto-haloes and their bound- 
aries in cosmological initial conditions (Porciani, Dekel & 
Hoffman, in preparation). 

The collapse proceeds in two succesive stages, corre- 
sponding to the quasi-linear and the non-linear regimes. The 
first stage is characterized by a two dimensional collapse in 
the plane defined by the major and intermediate inertia axes. 
This leads to the formation of a transient clumpy elongated 
structure aligned with the large scale filament within which 
the halo resides. In the fully non-linear regime the proto- 
halo experiences a one-dimensional collapse and a series of 
merger of the sub-clumps along the filament, leading to a 
quasi-spherical object in virial equilibrium. 

The proto-haloes can be characterized by a variety of 
other statistics, which may distinguish them from over- 
density patches about random density peaks. In terms of 
shape, most proto-haloes have significant ellipticities and 
most of them are prolate, but this does not distinguish them 
very clearly from over-density patches about peaks. In con- 
trast with the strong alignment of the inertia and deforma- 
tion tensors, there are only weak correlations between the 
triaxiality of the inertia tensor and that of the deformation 
tensor, as expressed by ratios of eigenvalues. 

The vast majority of proto-haloes are initially collaps- 
ing along two principal directions (except when the self grav- 
ity is also considered, inducing a collapse in all directions). 
This is in clear contrast with the behaviour of density peaks, 
which have similar probabilities to collapse along three or 
two axes. Our finding is consistent with proto-haloes being 
embedded in large-scale filaments surrounded by voids. The 
smoothed central densities in proto-haloes typically repre- 
sent la positive perturbations, but with smaller dispersion 
than for general peaks. This is partly because the proto-halo 
centers are only weakly correlated with the density peaks, 
and also because some of the high peaks lead to clumps that 
later merge into the big haloes we identify today and include 



in our sample. Despite the significant asphericity, the proto- 
halo densities agree on average with the predictions of the 
spherical collapse model, but the scatter is large. For exam- 
ple, there are proto-haloes with S <C I (with S the density 
contrast smoothed on the halo scale and linearly estrapo- 
lated to today), whose collapse is boosted by strong shear 
of large compression and dilation factors. 

In Paper I, we evaluated the performance of linear TTT 
in predicting the spin of galactic haloes. We found that, for 
a given proto-halo at the initial conditions, TTT provides 
a successful order-of-magnitude estimate of the final halo 
spin amplitude. The TTT prediction matches on average 
the spin amplitude of today's virialized haloes if linear TTT 
growth is assumed until about to/3. This makes TTT use- 
ful for studying certain aspects of galaxy formation, such as 
the origin of a universal spin profile in haloes (Bullock et al. 
2001b; Dekel et al. 200f), but only at the level of average 
properties, because the random error is comparable with the 
signal itself. We also found in Paper I that non-linear evo- 
lution causes significant variations in spin direction, which 
limit the accuracy of the TTT predictions to a mean er- 
ror of ~ 50°. Furthermore, spatial spin-spin correlations on 
scales > l/i _1 Mpc are strongly weakened by non-linear ef- 
fects. This limits the usefulness of TTT in predicting intrin- 
sic galaxy alignments in the context of weak gravitational 
lensing (Catelan et al. 200f; Crittenden et al. 2001). This 
situation may improve (as in the case of the spin-shear cor- 
relation addressed in the current paper) if the orientations 
of today's discs were determined by the halo spins at a very 
high redshift, which are better modelled by TTT. On the 
other hand, we know this only for the haloes that we se- 
lected at z = 0, and we don't know whether this is true for 
the subhaloes which host disc formation at higher redshift. 

In our studies of the cross-talk between the different 
components of TTT, we have also realized another surprise 
that leads to a revision in the standard scaling relation of 
TTT (White 1984). It is demonstrated in Paper III that the 
off-diagonal, tidal terms of the deformation tensor, which 
drive the torque, are only weakly correlated with the diag- 
onal terms, which determine the overdensity at the proto- 
halo centre. The latter enters the TTT scaling relation via 
the expected collapse time of the proto-halo, and the lack 
of correlation with the torque leads to a modification in the 
scaling relation. The revised scaling relation can be applied 
shell by shell in order to explain the origin of the universal 
angular-momentum profile of haloes (Bullock et al. 2001b; 
Dekel et al. 2001). It can also be very useful in incorporat- 
ing spin in semi-analytic models of galaxy formation (Mailer, 
Dekel & Somerville 2002). 
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